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The stability of the de Sitter era of cosmic expansion in spatially curved homogeneous 
isotropic universes is studied. The source of the gravitational field is an imperfect fluid 
such that the parameters that characterize it may change with time. In this way we 
extend our previous analysis for spatially-flat spaces as well as the work of Barrow. 

1. Introduction 

Inflationary cosmic expansions driven by a dissipative fluid has attracted some at- 
tention in the past as a mechanism able to solve the problems usually attached to 
the standard hot big-bang scenario (i.e. flatness, horizons, and monopoles abun- 
dance) different from those based on one (or more) scalar fields. The main dynamic 
effect of either the scalar field(s) or the dissipative fluid is to produce a state of 
cosmic accelerated expansion through the violation of the strong energy condition. 
If this accelerated state is sufficiently long, and if the Universe is able to exit it, 
the mentioned problems may find a satisfactory solution. Very often this analysis 
has been restricted to spatially-fiat Friedmann-Lemaitre-Robertson- Walker (FLRW 
for short) space-times, if only because the dependence of the scale factor on the 
spatial curvature k becomes negligible shortly after infiation has set in. However, 
one should be mindful that this curvature may have some impact on the stability 
of the inflationary solutions. Given the potential importance of such a scenario, it 
is advisable to stablish its generality with respect to the effect of spatial curvature. 
In El we studied the stability of inflationary expansions caused by a dissipative fluid 
governed by a causal transport equation (i.e. one that respects both relativistic 
causality and hydrodynamic stability of the fluid i, i), assuming fc = 0. Before 
going any further, a word of caution seems in order. The stability of the dissipative 
fluid, in the sense this expression is used here, depends on the physical .properties of 
the fluid under consideration, i.e. on its state and transport equations □, i. Accord- 
ingly it should not be confused with the stability of the cosmic expansion -either 
inflationary or not. The stability of inflationary expansions driven Jpy^dissipative 
fluids in curved FLRW space-times was briefly discussed by Barrow Q, However, 
his study was confined to fluids governed by transport equations of Eckart's type 
which, as is well-known, entails the aforesaid problems of acausality and instability. 
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The aim of this paper is to investigate the combined effects of viscosity and 
curvature extending the analysis of EJ to spatially curved FLRW space-times and to 
study scenarios not considered previously. This bears some interest because obser- 
vationally it is doubtful whether the spatial curvature of the Universe is positive, 
null or negative. On the one hand the average matter density is well below the 
critical value (~ 0.3 in critical units), but on the other hand the Universe seems to 
have entered an accelerated phase nowadays, which strongly hints to a cosmological 
constant of about 0.7 in the same units -however its value is limited by gravitational 
lensingi. Further, although the location of the first acoustic peak in the CBRgpec- 
trum seems to suggest a flat Universe, its exact position is still uncertain E£l, liil. In 
extending the results of we also extend the work of Barrow's since the transport 
equation we consider is of causal type and therefore more general than the one used 
by him -for a bjdfiOlkistLationxif the impact of causal transport equations in FLRW 
cosmology see Il3 Ej llJ lla and 

To analyse the asymptotic stability of the inflationary expansions we will resort 
again to the second method of LyapunovEO, succintly stated in the appendix of El. 
In addition, some of the scenarios considered in this paper will be similar to those 
of u, whereby we shall avoid unnecessary repetition of details. As it turns out, the 
set of de Sitter solutions is drastically restricted because of a constraint that links 
some of the fluid parameters with the rate of cosmic expansion. Our choice of units 
is c = SttG = /cb = 1. 

2. Asymptotic Stability 

We consider a FLRW universe filled with a imnerfectfluid whose dissipative bulk 
pressure obeys the causal transport equation H, llil, ta 

TT + Til = -'iQH (1) 

where r(> 0) is the relaxation time associated at the dissipative process, C(> 0) the 
phenomenological coefficient of bulk viscosity, H = a{t)/a{t) denotes the Hubble 
function, and a{t) the scale factor of the FLRW metric. The corresponding Einstein 
field equations are 

H"^ = \p-^ (A = 0,±1) (2) 
3(ij + ij2)^_i[p + 3(p + ^)] (3) 

with p and p the energy density and hydrostatic pressure, respectively. These latter 
two quantities are assumed to be linked by the equation of state p — (7 — where 
7 denotes the adiabatic index of the cosmic fluid. This, in general, depends on time 
-though oftenly, for simplicity, it is assumed constant. 
The above set of equations yields 
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[(1-|7) (2i/ - T-i) + I7] - 



(4) 



and the latter can be recast as 



dt 



D{H,H). 



(5) 



Here, the left hand side is the time derivative of a Lyapunov function (see appendix 
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As shown in figure 1 the potential (^) has two extrema, a minimum at and a 
maximum at . Their asymptotic behaviours for an expanding evolution at late 
time are 



H+ = 



c 



7 + r7 



0{a~ 



and 



H- 



7 T 



(8) 



(9) 



The solution (H) is unstable, while the stability of the solution is determined by 
the sign of the leading term of _D in a neighbourhood of the point (iJo,0) of the 
phase plane {H, H) , where i7o is the leading term of i7+ 



D^- {2,^H + r-i) + (^i/ - I (37i? + r-^) + 0{a-'). (10) 

Here we have assumed that i?+ is a quasi-de Sitter solution. A sufficient condition 
for this solution to be asymptotically stable is that 87^ + r^^ increases for large 
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time. This is equiyalcnt to the condition that <^/t increases with time. Then using 
the relationship lij 

^ = v^^p (11) 
r 

where v is speed of the dissipative signal associated to tt, this condition is also 
equivalent that ^v^ is an increasing function of time. 

It can be seen that Hq is an exact de Sitter solution when the constraint 

7 + [r-i - 2i7o] [t-^ - 2i?o] (12) 

is satisfied. It gives for 7 7^ 2/3 



° 2(37-2)t ^ ' 



and 



H, - (14) 

when 7 = 2/3. Also the static solution H ^ exists for A: = 1, when the constraint 

2 , , 

T7 + 7=- (15) 

holds. To analyse its stability we linearize equation (0) about it, 



jr 1 
^ " 2 



H ^0. (16) 



The general solution of ( p^ ) grows for any decreasing adiabatic index showing that 
the static solution is unstable and can be indentified with the limiting case of the 
unstable solution This result could be different if a non-standard increasing 
adiabatic index in the limit t — > cxd were admitted. However, by fine tuning the 
initial condition a one-parameter family of solutions that approaches a static space- 
time at large times can be obtained. Besides, there exists a one-parameter family 
of solutions that starts from a static spacetime in the far past and evolves towards 
a stable de Sitter solution. 



3. Inflationary Solutions 

Accelerated expansion requires a sufficiently negative pressure, hence we will con- 
sider evolutions that approach asymptotically to a negatively constant viscous pres- 
sure. To this end we will investigate three models, two of them have the viscous 
pressure coefficient proportional to Hubble time, and in the last one it is determined 
by causality and stability of scalar perturbations. These three cases are character- 
ized by the specific instances: (i) when 7 = and C cx H~^, (ii) when r and C vary 
as H~^, and (iii) when oc H and C = v^^pr. We will see them in turn. 
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3.1. Constant adiabatic coefficient model 

When C = constant (0) can be rewritten as 



dt 



= 



accordingly, the bulk viscosity coefficient must be given by 

c 

with c a positive constant, and 



2c 



(17) 



(18) 



(19) 



is a first integral of (|17|). Assuming that the adiabatic index does not vary with 
time, this equation can be easily reduced to an equivalent mechanical system by the 
substitution a — s^/'^^ (in the appendix we obtain its parametric general solution), 
with first integral 



^ 97 fc 2(1-2/37) 

2+8 



37c 2 



3jE 



(20) 



where E plays the role of a mechanical energy. This expression shows that the 
curvature and viscosity terms behave as "conservative forces" . In particular vis- 
cosity provides an effective spring with negative elastic constant, and it determines 
the late time behaviour of this model. Next we will integrate thjs system for some 
values of 7 that yield simple solutions and were considered in □ to compare with 
the results obtained for fluids governed by Eckart's transport equation: 
I) When 7 = 1/3 wc obtain 



a{t) = ci + C2 exp 



provided the constraint 8CC2C3 — k — 2cc\ 
perfect fluid solutions are recovered 



2ct ) + C3 exp i^-Vi'cVj (21) 
holds. In the limit c — > the usual 



a{t) = ait^ + a2t + 03 (22) 

provided the constraint 4aia3 ~ ~ k — Q is satisfied. This family of solutions 
includes several interesting possibilities. Solutions with and without initial singu- 
larity, bouncing solutions and solutions with a finite time-span, as well as the Milne 
solution when fc = — 1 and ai = 0. 

II) When 7 = 2/3 there are three families of solutions depending on the sign of 

E: 
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1. For E>0 

a{t) = W — sinh ( VcAi) (23) 

i.e. the scale factor initially grows as \j2Etd independently of c. 

2. For < 

/ 2^ 

a{t) ^ J cosh (VcAt) (24) 

this family of non-singular solutions bounce at At = when the scale factor 
attains its minimum value = ^ j c. 

3. For E = ^ 

a{t) = y^exp {y/^At) . (25) 

The three families of solutions exist for all fc = or ±1, and for large time they 
approach the inflationary de Sitter scenario a{t) w exp ^/cAt. These results extend 
those of §4 in U obtained in the framework of the non-causal Eckart's theory, and 
where the viscosity was assumed to follow a power-law dependence upon the density. 
Ill) When 7 = 4/3 (i.e. a radiation-dominated universe) one has 



fc + fcexp (-\/2^i) ^h)^^^ '^^^^ {V^'ct 



(26) 



where h is an integration constant. This familiy of solutions comprises three types 
of behaviours depending on the value of E. When E > singular solutions occur 
for all k. However, when — l/4c<i?<0a singularity arises only if fc < 0. Near 



the singularity the behaviour of the solution (26) is radiation-like. When E < —c 
there is no singularity. 

Again the large-time behaviour is qualitatively the same as that found with 
Eckart's theory. Otherwise the evolution behaves differently. 

IV) When 7 > 2/3 it follows, using the results of the appendix, a family of 
solutions that near the singularity evolves as a « t^/'^'^ , and at late time it exhibits 
a de Sitter expansion a w exp ■\/2c/37t, irrespective of the spatial curvature. 

3.2. Constant interactions number per expansion time model 

The expanding universe defines a natural time-scale - the expansion time . 
Any particle species will remain in thermal equilibrium with the cosmic fluid so 
long as the interaction rate is high enough to allow rapid adjustment to the falling 
temperature. The relaxation time t is the characteristic colisional time of hydro- 
dynamical procesess occurring after the quantum era. Then a necessary condition 
for maintaining thermal equilibrium is 

T < R-^ (27) 

Now r is determined by 



6 



1 

nav 



(28) 



where a is the interaction cross-section, n is the number density of the target 
particles with which the given species is interacting, and v is the mean relative speed 
of interacting particles. So we have that — (tH) is the number of interactions 
in an expansion time. Now we will consider the case when is a constant larger 
than one. 

To find the solutions of equation (^ that satisfy this property we insert the 
ansatz 



7= 3(1 + ^) 



(29) 



along with (|lj) in the equation (||). A set of solutions can be found for a variable 
adiabatic index. Expression splits in two equations. 



H +{2 + iy)HH + vH^ - ucH = 
plus a linear equation in e. 



2H H 



i/2 



provided r ^ = vH with 1/ a constant. Solving (|3l]), we get 
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1 



(30) 



(31) 



(32) 



where 6 is a free parameter. To make explicit the t dependence of the adiabatic 
index we need solutions of ( ^0[ ) that satisfy < 7 < 2. Equation (^0|) can be 
transformed into a second order linear differential equation whose general solution, 
obtained in Ei, reads 

H^^c+% + ^ (33) 



where ci and C2 are arbitrary integration constants. There it was shown that simple 
explicit solutions of (^0|) can be obtained when v — 1 and when = 4. In this case 
the solutions are 



i{t) = ci + C2 exp (Vci) + C3 exp (-%/ci) 



and 



(34) 
(35) 



a{t) = [ci + C2 exp {2^/ct) + C3 exp [-2^/ct)] 

respectively. For t ^ c^^/^, both ( |3^ ) and (|35| ) describe inflationary expansions 
regardless the initial conditions. 



3.3. Stable causal sound perturbations model 
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In this subsection we will investigate the consequences imposed by causality and 
stability on scalar longitudinal sound perturbations. It can be seen that the re- 
lationship between the viscosity coefficient and the dissipative contribution to the 
sound speed is C = w^7pr. 

In this case equation (|^) transforms into 

h" + (37 + iy)h' + 37 {ly - 3v^) h = (36) 

where the variable h = + k/o? is proportional to the energy density of the fluid 
and the prime indicates derivative with respect to 77 = In a. Except for v — 3w^, 
the general solution of this equation leads asymptotically to power-law behaviors 
When 1/ = the solution exhibit de Sitter expansion at late time. Exact 



solutions of m6f) are 



a{t) = sinh (VbAt^ , (fc -1) (37) 

a{t) = cosh (VbAt^ , (fc = 1) (38) 

where b is an integration constant. These solutions are attractors and therefore of 
singular importance as they indicate the leading behaviour for large cosmological 



times. It is remarkable that the expressions (|3^) and (38) do not explicitly contain 
any quantity directly associated to viscosity (or particle creation) even though it 
is precisely this effect that drives the exponential expansion. In this case the same 
effect generates the translational invariance of the energy density. 

4. Additional Inflationary Scenarios 
4.1. Variable r 



When either 7 or C is a constant the solution (13) is asymptotically stable provided 
T is a decreasing function. When r is a function that first decreases and in a second 
stage increases, there is a first period of exponential inflation followed by a graceful 
exit. On the other hand, when the relationship (|l^) holds, r presents a minimum 
provided has a maximum. In this case the Universe enters an inflationary stage 
and afterwards exits it. This parallels the cosmic scenario of §2.2 in lil. There the 
cosmic fluid was modelled as a mixture of radiation and heavy particles that decayed 
at a very high rate into (more_stable) lighter particles (less massive modes), with 
high or moderate multiplicity t3. 

4.2. Variable t and 7 

Another model of dissipatively driven inflation arises when is a constant (or nearly 
a constant) and 7 varies as 



7(0 4 



c 



(39) 
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with c a constant. From we see that D is negative-definite when f < 0, and 
accordingly the de Sitter expansion results asymptotically stable. Likewise 7(t) 
increases when c > 0. 

By choosing c = T(ii)/2 and T{ti)/T(t2) = 2 one follows a model in which 
initially all the energy is in the form of non-relativistic particles, 7 = 1, and it is 
gradually transferred to relativistic ones, so that finally 7 = 4/3. This corresponds 
to the dissipative process of decay of dust particles into radiation. Additional dis- 
sipation may arise from the interaction matter-radiation. Once the heavy particles 
have decayed, the Universe exits inflation. A similar scenario in flat space was re- 
ported in El, (see equations (17) & (18) there) but in that work r was considered 
constant instead. 

Assuming that the dust particles are primeval mini-black holes of rest mass m 
and that the thermodynamic properties of the mixture of black holes and relativis- 
tic particles correspond to a Boltzmann gas of vanishing chemical potential, the 
adiabatic index reads 

^(^) = 1 + ^^^^^^ ^ ^ (40) 

where are modified Bessel functions of the second kind, and z = m/T the 
dimensionless inverse temperature. By equating the right hand sides of ( ^9|) and 
( ^0|) we can describe the continuous process of decay of mini-black holes from t — ti, 
when the black hole energy density dominates the Universe, until t2 when the mini- 
black holes have completely evaporated away and the Universe becomes radiation- 
dominated. Here it is understood that all the black holes have the same mass and 
therefore the same temperature, and that this one equals the temperature of the 
massless component of the cosmic fluid aLthe beginning of the evaporation. The 
temperature behaves as in the A: = case EJ 

z[ ^ 12K2{z) + 3zKi{z)- Bz^ 

z ^ 12K2iz) + 5zKi{z) + z^Ko{z) ^ ' 

where B = 9H(^/Aq, with a constant and ' = d/Hdt. From this equation follows 
that z' is negative for large z, which agrees with the warm inflationary scenario of 
above. 

The time dependence of this temperature near ti and t2 for a generic T{t) can 
be made explicit by expanding ( |39| ) about t = ti and (^) for z — > oo. Thus one 
follows T cx t — ti, while in the opposite limit Ji.e. t t2 and z — > 0) it yields 
T (X {t2— t)^^/^. This reproduces the results of EJ but this time with a generic rit). 
Again the production of relativistic particles at the final stage of the black holes 
evaporation is accompanied by a huge increase of the tempemture of the cosmic 
fluid -which also agrees with a previous study of this process E3. 

The entropy productianjier unit volume in the radiation fluid is given by the 
well-known expression El, til, 5 = n^/{(T). It has the limiting behaviours 

t^U (42) 
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where 



and 



viC, 



t ~* t2 



1 1 



I 2r(t) 



3r(ti)C(2 + -r(ti)C) + 



1 - (1 + 3T(ii)C) v/l + 6r(ti)C + 1 



+ 



+ ^) (l + 3r(ii)C - v/l + 6r(ii)c) + 



a{t) — ao exp ■ 



v/l + 6Cr(ti)~l 

2T(tl) 



(43) 



(44) 



(45) 



The entropy production rate in the radiation fluid happens to be very high at the 
beginning of the evaporation, but decreases sharply about the end of this process. 
As in B the net rate of radiation particle production per mini-black hole and unit of 
volume varies roughly as (p + p)^^, where in this case p and p refer to the radiation 
fluid only B. 



5. General case 

In this section we consider the de Sitter solution (|3|) and assume that 7, r and C 
vary with time. Again this more general situation may occur during the decay of 
massive particles into lighter ones and also during the decay of four-dimensional 
fundamental strings into massive and massless particles -admittedly this second 
possibility is more speculative. In this case the algebraic relationship holds 



2Hn 



3^ 6H§T{t) 



(46) 



By adequately choosing the behaviour of C(<) and r(t) the decay into radiation, 
reaching relativistic gas state when 7 = 4/3, can be described. After the decay a 
condensation phase back into non-relativistic matter may occur. It ends when 7 
returns to 1. A scenario corrmatible with the latter phase is the quantum tunneling 
of radiation into black holes E3. This may arise very naturallvjiecause of the insta- 
bility of the hot radiation against spontaneous condensation c3. (This is altogether 
different from the whole disappearance of the radiation by black hole accretion). 
During this period both the viscosity coefflcient and the relaxation time may be 
chosen as monotonic decreasing functions, provided C/r grows with time. 

Again we may interpret this behaviour in terms of a two-fluid model, where the 
viscosity coefficient arises because of the particle production process from the decay 
of massive non-relativistic particles into light ones. Shortly after the beginning of 
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the decay the particle production rate is large and the energy density of the fluid 
becomes dominated by the light component. Later on, as the decay rate slows down, 
the effect of adiabatic dilution by the fast exponential expansion of the Universe 
turns out to be more important. Accordingly the non-relativistic energy density 
takes over again, since it goes down as a~^, while the relativistic energy density 
goes down at the faster rate of a~^. 

The adiabatic index can be used to estimate the dissipative contribution to 
the speed of sound in (1l]l) 



= ^^"^ (47) 



In the limit t —> oo, is a monotonic function of C and t within the region 2/3 < 7. 
It vanishes when C = a-nd reaches a maximum value v"^ = 2/3 when C/^o ^ 1- 

6. Concluding Remarks 

We have applied the second method of Lyapunov to analyse the stability of cosmic 
inflationary expansions, driven by a dissipative fluid governed by a transport equa- 
tion that allows for relaxation (i.e. of causal type), for non spatially flat FLRW 
metrics. The parameters characterizing the fluid (adiabatic index, viscosity coeffi- 
cient and relaxation time) may vary as the expansion proceeds. This is interesting 
as repeteadly stressed, a dissipative fluid can phcnomenologically mimic a perfect 
one where particle production takes place i, i, i, 0, either from the quantum 
vacuum, or by the decay of pre-existing heavy particles, or by the decay of massive 
modes of fundamental strings into massless modes. 

Essentially the curvature does not modify the results found in our previous paper 
EJ concerning the stability of the de Sitter solution. This result shows the insensibility 
of the viscosity-driven inflationary scenario with respect to initial conditions even 
when spatial curvature is present. However, it severely restrics the set of de Sitter 



solutions encompassed by this model because a constraint (12), linking the adiabatic 
index with the relaxation time, must be satisfied. Nevertheless, since the de Sitter 
solution is asymptotically stable for a wide set of behaviours of the fluid parameters, 
inflation appears rather natural. The static solution, = 0, is found to be unstable; 
however this result may be altered if a non-standard increasing adiabatic index could 
enter the jjilay. This work generalizes p,wn previous study in spatially-flat FLRW 
universes □ as well as those of Barrow □, El. Further, new solutions describing the 
effects of dissipation have also been found. They either go over a stable de Sitter 
scenario from an initial power-law singularity or asymptotically approach a static 
universe if a fine tuning of the initial condition is made. 

Lastly, an expression for the dissipative contribution to the sound speed have 



been obtained (47). In general the latter goes down as the Universe expands, 
something rather natural as the impact of disipative effects is thought to diminish 
with expansion. 
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Appendix 

Equation ( p^ with the change of variable a = s^/^^ becomes 
s + F{s)^Q^ F(.) = 1^ - l) ^.1-4/3-^ 
Then, by the nonlocal transformation 

../fW*, ,^/f(„)., (49) 
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equation ( pq ) turns into a linear second order difFcrcntial equation with constant 
coefficient 

z" + 1 = (50) 

where the prime indicates derivative with respect to rj. Solving it the general so- 
lution of ( |4^ ) in parametric form follows. This can be achieved by inserting this 
solution in the second equation of (|49| ) 

-^+c,, + c. = ^fc.^(-^/3^)-^.^ (51) 

The transformation law between rj and t follows from solving (|^) for s and inserting 
the resulting expression in the second equation of ( ^9|) 



Figure Caption 



Figure 1 Graph of the potential V{H), given by equation (^), for each value of the 
spatial curvature k. 
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